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Abstract. The masses of charmonium s and /^-states, pseudoscalar and vector decay constants, lep- 
tonic, hadronic as well as radiative decay widths for charmonia have been computed in the frame- 
work of extended harmonic confinement model without any additional parameters. The outcome in 
comparison with other contemporary theoretical and experimental results is presented. 
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The experimental investigations in charmed quark sector has opened up a vast data 
bank for the hadronic decays into various channels [[ll, SlU]. A hadronic transition or 
decay between quarkonium levels is understood as a two-step process in which gluons 
first are emitted from the heavy quarks annihilating each other and then recombine 
into light hadrons. Perturbative QCD is not directly applicable because the energy 
available to the light hadrons is small and the emitted gluons are soft. Nevertheless, 
the final quarkonium state is small compared to the system of light hadrons and moves 
nonrelativistically in the rest frame of the decaying quarkonium state. 

Decay rates of heavy quarkonium states in the non-relativistic limit with some QCD 
radiative correction have been studied through various phenomenological models sub- 
jected to substantial relativistic corrections [4]. In the heavy flavour sector, the nonrela- 
tivistic expressions in accordance with the factorizable assumption for the long-distance 
part and the short-distance part of the decay amplitudes, are proportional to the square 
of the radial wave function at the origin. We employed there a relativistic harmonic 
confinement model with scalar plus vector potential 150 . 



HADRONIC MASSES FROM ERHM 

The mass of a hadron having p number of quarks in ERHM can be obtained as Isj], 
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where the first sum is the total confined energies of the constituting quarks of the hadron, 
the second sum corresponds to the residual colour coulomb interaction energy between 
the confined quarks and the third sum is due to spin dependent terms. 



The intrinsic energies of the quarks in charmonium system is given by 
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The coulombic part of the energy is computed using the residual coulomb potential 
given by [6], 
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here (0„ represents the colour dielectric "coefficient" |@]. (0„ is found to be to be state 
dependent fsl], so as to get consistent coulombic contribution to the excited states of 
the hadrons. It is a measure of the confinement strength through the nonperturbative 
contributions to the confinement scale at the respective threshold energy of the given 
flavour production. 

We construct the wave functions for charmonia by retaining the nature of single parti- 
cle wave function but with a two particle size parameter Q.M{qiqj) instead of Q.N{q) I0]- 
Now coulomb energy is computed perturbatively using the confinement basis with two 
particle size parameter defined above for different states as, £N{qiqj)coui = {N\Vcoui\N) . 
The fitted parameters to obtain experimental ground state mass are rtic = 1428 MeV, 
k = 0. 19252 and the confinement parameter A = 2166 MeV^/^. 

The confined gluon propagator derived in the current confinement scheme (CCM) 
for gluons has been employed in the derivation of one gluon exchange potential using 
RHM wave functions for quarks under the Breit formalism ISQ. This COGEP has central, 
tensor and spin-orbit terms in addition to many momentum dependent terms. 

From the centre of weight masses, the pseudoscalar and vector mesonic masses are 
computed by incorporating the residual two body chromomagnetic interaction through 
the spin-dependent term of the COGEP perturbatively as, 

eUqiqj)s.D. = {NJ\Vsd\NJ) (4) 

where |A^7) is the given hadronic state. For mesons |A^7) becomes the Iqiqj) states. We 
consider the two body spin-hyperfine interaction of the residual (effective) confined one 
gluon exchange potential (COGEP) as 

The spin-orbit interaction of the residual (effective) confined one gluon exchange poten- 
tial (COGEP) is given as LSJ 
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In the case of mesonic states, L and S correspond to the orbital and spin operators of 
the relative co-ordinates. Where DQ^Vij) and D[{rij) appeared in Eqns ([5]) and © are 
derivatives of the confined gluon propagators of CCM [[SO. 



TABLE 1. Masses of cc in comparison with PDG and other theoretical models 
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DECAY PROPERTIES OF MESONS 

Without any additional free parameters, using the Van-Roy en- Weisskopf formula for 
computing leptonic decay widths [|9|] without radiative correction term, we have com- 
puted leptonic decay widths and are tabulated in Table [2] alongwith other theoretical as 
well as experimental values. 

The Van Royen - Weisskopf formula used for the meson decay constants is obtained 
in the two-component spinor limit [[loll . fp and fy are related to the ground state radial 
wave function 7? 15(0) at the origin, by the Van Royen - Weisskopf formula modified 
for the colour as, f^^y = 3\RisiO)\'^ / TlMp/y where Mp/y is the ground state mass of 
the pseudoscalar/vector meson. In the present computations, fp = 410.34 MeV which is 
closer to the experimental value of 410 ± 15 MeV, while the other results are 496.42 



!D, 357 UM, 349 1120, 498.90 [13], 549.61 [14] and 663.06 [15] (all values in MeV). 
We employ radial wave functions to compute the hadronic as well as radiative decay 
widths of the vector and pseudoscalar mesons of cc system based on the treatment of 



perturbative QCD as 
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Here, symbols have their usual meanings. The computed decay widths in comparison 
with other theoretical predictions are tabulated in Table [3l 



TABLE 2. Leptonic decay widths (in keV) of cc{n^Si ) 
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TABLE 3. Tjc ^ YY (in ke V) with correction terms and light hadron decay (in MeV) 
of cc meson 
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